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Abstract 

Deep-inelastic positron-proton scattering events at low photon virtuality Q 2 with a for- 
ward jet, produced at small angles with respect to the proton beam, are measured with the 
HI detector at HERA. A subsample of events with an additional jet in the central region 
is also studied. For both samples differential cross sections and normalised distributions 
are measured as a function of the azimuthal angle difference, Acj), between the forward 
jet and the scattered positron. The sensitivity to QCD evolution mechanisms is tested by 
comparing the data to predictions of Monte Carlo generators based on different evolution 
approaches as well as to next-to-leading order calculations. 
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1 Introduction 



Measurements of the hadronic final state in deep-inelastic lepton-proton scattering (DIS) test 
Quantum Chromodynamics (QCD), the theory of the strong force. At moderate negative four- 
momentum transfers squared Q 2 of a few GeV 2 , the HERA ep collider has extended the avail- 
able kinematic range for deep-inelastic scattering to regions of small Bjorken-x ~ 10 -4 . This 
is the region of high parton densities in the proton, dominated by gluons and sea quarks. At the 
large 7*p centre-of-mass energy available at small x, a transition is expected from parton cas- 
cades ordered in transverse momentum, described by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equations [1], to cascades unordered in transverse momentum, de- 
scribed by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach . 




Figure 1: Generic diagram for deep-inelastic ep scattering at small x. A gluon cascade evolves 
between the quark box, attached to the virtual photon, and the proton. The gluon longitudinal 
momentum fractions and transverse momenta are labeled xi and k Ti , respectively. 

A generic diagram for parton evolution in a DIS process at low x, in which a gluon from the 
proton induces a QCD cascade before an interaction with the virtual photon, is shown in figure[[] 
In the DGLAP approximation the struck quark originates from a parton cascade ordered in 
virtualities of the propagator partons. At low x this implies a strong ordering in transverse 
momentum, k T , of the emitted partons, measured with respect to the proton direction. In the 
BFKL approach there is no ordering in kr of the partons along the ladder. Compared to the 
DGLAP scheme more gluons with sizable transverse momentum are emitted near the proton 
direction. For this reason energetic jets of high transverse momentum produced close to the 
proton direction in the laboratory frame, referred to as the forward region, are considered to 
be especially sensitive to QCD dynamics at low x 0. Forward jet production was measured 
previously by the HI and ZEUS collaborations. In these measurements as well as in the present 
one, the requirements on the forward jet and the phase space were chosen in such a way that 
the standard DGLAP evolution is suppressed and the effects of BFKL dynamics are enhanced. 
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Preference for models which employ QCD evolution non-ordered in transverse momentum was 
observed PHIDI. 

One of the observables suggested to be sensitive to BFKL dynamics ifTTI is the azimuthal 
angle difference, A0, between the forward jet and the scattered electron, defined in the labo- 
ratory frame. In the Quark Parton Model (QPM) process e + q — > e + q the simple two-body 
kinematics constrains the scattered electron and the jet to be produced back-to-back, and thus 
predicts at the parton level A</> = it. Hadronisation effects induce some smearing to this parton 
level prediction. Inclusion of higher order processes partially decorrelates the jet from the elec- 
tron. As a consequence, for evolution schemes without ordering in transverse momentum, the 
decorrelation is expected to increase with electron-jet rapidity distance, Y, since the phase space 
for additional parton emissions increases. The calculations employing the BFKL approach to 
the next-to-leading order accuracy (NLO BFKL), indeed predict an increase of the azimuthal 
angle decorrelation with the electron-jet rapidity distance lfl2ll . 

This paper presents a study of low x DIS interactions in which high transverse momentum 
jets are produced in the forward region. The forward jet cross sections and normalised distri- 
butions are measured as a function of the azimuthal angle difference Acf) in three bins of the 
rapidity separation Y between the positron and the forward jet. The forward jet cross section as 
a function of Y is also measured. Moreover, the measurements of the azimuthal correlations in 
A0 are performed using a subsample defined by a requirement of an additional central jet. In 
comparison with the forward jet sample, this subsample is expected to contain a higher fraction 
of forward jets from additional gluon emissions. 

The data set used for the analysis was collected with the HI detector in the year 2000, 
when positrons and protons collided with energies of 27.6 GeV and 920 GeV, respectively, 
corresponding to a centre-of-mass energy of y/s = 319 GeV. The integrated luminosity of the 
data set is 38.2 pb _1 , which is about fourteen times larger than in the previous measurement of 
the azimuthal decorrelation of forward jets Q. 



2 QCD Calculations 

The measurements presented here are compared with predictions of Monte Carlo (MC) pro- 
grams and perturbative QCD calculations at next-to-leading order (NLO). The MC programs 
use first-order QCD matrix elements and model higher order terms by parton showers in the 
leading logarithm approximation or by quasi-classical gluon radiation from colour dipoles. 
Three MC event generators, which adopt different QCD based approaches to model the par- 
ton cascade, are used. 

• RAPGAP [13J matches first order QCD matrix elements to DGLAP based leading-log 
parton showers with kx ordering. The factorisation and renormalisation scales are set to 
fif = fi r = \/Q 2 + p\, where is the transverse momentum of the two outgoing hard 
partons in the centre-of-mass of the hard subsystem. Predictions of RAPGAP are labeled 
DGLAP in the figures. 
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• DJANGOH HH with ARIADNE includes an implementation of the Colour Dipole Model 
(CDM) [15 1, which has as its basic construct a colour dipole formed by the struck quark 
and the proton remnant. Subsequent parton emissions originate from a chain of indepen- 
dently radiating dipoles formed by the emitted gluons. In this approach the transverse 
momenta of emitted gluons perform a random walk such that CDM provides a BFKL- 
like approach. The leading order partonic final state is corrected to exactly reproduce the 
0(as) matrix elements. The simulation of DJANGOH/ ARIADNE uses a set of colour 
dipole parameters tuned to describe measurements of the hadronic final state in DIS at 
HERA QU. The DJANGOH/ARIADNE predictions are referred to as CDM in the fol- 
lowing. 

• CASCADE [HI implements the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolu- 
tion [fT8l which aims to unify the DGLAP and BFKL approaches. It introduces angu- 
lar ordering of emissions to implement gluon coherence effects, and thus in the high 
energy limit the CCFM evolution equation is almost equivalent to the BFKL approach, 
while reproducing the DGLAP equations for large x and high Q 2 . CASCADE uses off- 
shell leading order QCD matrix elements, supplemented with gluon emissions based on 
the CCFM evolution equation, requiring an unintegrated gluon density function (uPDF) 
which takes the transverse momenta of the propagators into account. In this paper two 
different uPDF sets are used: set AO |fP9l with only singular terms of the gluon splitting 
function and J2003-set 2 [|20l including also non-singular terms, labeled set 2 in the fig- 
ures. These parameterisations for the unintegrated gluon density were obtained using the 
CCFM evolution equation to describe the structure function F 2 (x,Q 2 ) as measured by 
HI [EQ and ZEUS (221. Predictions of CASCADE are labeled CCFM in the figures. 

To perform the hadronisation step, all of the above models use the Lund string fragmenta- 
tion scheme, as implemented in JETSET [|2J| in the case of DJANGOH/ARIADNE and in 
PYTHIA for RAPGAP and CASCADE, using a tuning based on LEP e+e~ data E3. The 
RAPGAP and DJANGOH/ARIADNE predictions are calculated using the HERAPDF1.0 ||2f0 
set of parton distribution functions (PDF). 

The RAPGAP and DJANGOH/ARIADNE programs are interfaced with HERACLES B27H. 
which allows the simulation of QED-radiative effects. These MC models are used to simulate 
detector effects in order to determine the acceptance and efficiency for selected forward jet 
events in DIS. Generated events are passed through a GEANT ll28l based simulation of the HI 
apparatus, which takes into account the running conditions of the data taking. Simulated events 
are reconstructed and analysed using the same program chain as is used for the data. 

The measurements of azimuthal correlations are also compared to the fixed order NLO 
DGLAP predictions of NLOJET++ ||2~9ll . The NLOJET++ program is used here to calculate 
dijet production at parton level in DIS at NLO(a|) accuracy. It should be noted that the jet 
search is performed on partons in the Breit frame (see section 3.2), and therefore the events 
contain at least one jet in addition to the forward jet. The renormalisation and factorisation 

scales are defined for each event and are set to \i r = [if = ^/(-P| sc + Q 2 )/2, where P T ,sc 
is the transverse momentum of the forward jet or the average transverse momentum of the 
forward and central jet in the forward jet sample and in the sample with an additional central jet, 
respectively. The NLO calculations are performed using the CTEQ6.6 ll30l parameterisation of 
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the parton distributions in the proton. The NLOJET++ parton level cross sections are corrected 
for hadronisation effects using the RAPGAP model. The correction factors for hadronisation 
are estimated bin-by-bin by calculating the ratio between the cross section for jets reconstructed 
from stable hadrons (hadron level) and the parton level cross section. The correction factors 
for hadronisation are in the range from 0.90 to 1.08, increasing with rapidity distance Y. The 
uncertainty of the NLOJET++ predictions due to missing higher orders is estimated by applying 
a factor 2 or 1/2 to the renormalisation and factorisation scales simultaneously. 



3 Experimental Method 
3.1 HI detector 

A detailed description of the HI detector can be found elsewhere ll3"lT[3~3~1l . The components of 
the detector which are most relevant for this analysis are briefly described below. The origin 
of the HI coordinate system is the nominal ep interaction point. The direction of the proton 
beam defines the positive z-axis. Transverse momenta p T and polar angles 9 of all particles are 
defined with respect to this direction. The azimuthal angle <p defines the particle direction in the 
transverse plane. The pseudorapidity is given by r] = —In (tan 9/2). 

The ep interaction region is surrounded by the central tracking detector (CTD) consisting of 
two large concentric drift chambers, operated inside a 1.16 T solenoidal magnetic field. Charged 
particles are measured in the angular range 20° < 9 < 160° with a transverse momentum 
resolution of cr PT /p T w 0.005 ■ p T [GeV] © 0.015. Information from the CTD is used to trigger 
events, to locate the event vertex, and contributes to the reconstruction of the hadronic final 
state. 

A highly segmented liquid argon (LAr) calorimeter is used to measure the hadronic final 
state. It covers the range of the polar angle 4° < 9 < 154° and offers full azimuthal coverage. 
The LAr calorimeter consists of an electromagnetic section with lead absorbers and a hadronic 
section with steel absorbers. The total depth of both sections varies between 4.5 and 8 inter- 
action lengths in the region 4° < 9 < 128°, and between 20 and 30 radiation lengths in the 
region 4° < 9 < 154° increasing towards the forward direction. Test beam measurements of 
the LAr calorimeter modules showed an energy resolution of <Je/E 0.50/ \/ E[GeV] © 0.02 
for charged pions (3l and of a E /E w 0.12/ y/E[GeV] © 0.01 for electrons fell . 

A lead/scintillating fiber calorimeter (SpaCal) 11331 covers the region 153° < 9 < 177.5°. 
It has an electromagnetic and a hadronic section and is used to measure the scattered positron 
and the backward hadronic energy flow. The energy resolution, determined from test beam 
measurements 11361 . is cte/E ps 0.07/ ^/ E[GeV] © 0.01 for electrons. The precision of the 
measurement of the polar angle of the positron, improved using the backward drift chamber 
(BDC) situated in front of the SpaCal calorimeter, is 1 mrad. 

The luminosity determination is based on the measurement of the Bethe-Heitler process 
ep — > ep'j where the photon is detected in a calorimeter located at z = —103 m downstream of 
the interaction region in the positron beam direction. 
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3.2 Event selection 



DIS events are selected using triggers based on electromagnetic energy deposits in the SpaCal 
calorimeter and the presence of charged particle tracks in the central tracker. The trigger effi- 
ciency is determined using independently triggered data. For DIS events with a forward jet, the 
trigger efficiency lies between 60% and 80%, and for the topology 'forward and central jet' it is 
at the level of 80%. 

The data set is restricted in inelasticity y, photon virtuality Q 2 and x: 0.1 < y < 0.7, 
5 < Q 2 < 85 GeV 2 , 0.0001 < x < 0.004. In this analysis these variables are determined 
from measurements of the scattered positron energy and its polar angle, and from the incident 
positron beam energy. This phase space is chosen to ensure that the DIS kinematics are well 
determined and to reduce the background from photoproduction. 

The background from photoproduction and from events with large initial-state QED radia- 
tion is further reduced by requiring 35 < — p z<i ) < 70 GeV. Here Ei and p z ^ are the 
energy and longitudinal momentum of a particle i, respectively, and the sum extends over all de- 
tected particles in the event. Energy-momentum conservation requires that Ej (Ei —p z ,i) = 
where E® is the positron beam energy. Jets are identified from combined calorimeter and track 
objects ||37l using the k T cluster algorithm in the longitudinally invariant inclusive mode |[38l 
applied in the Breit frame. The reconstructed jets are then boosted to the laboratory frame. 

The measurements of forward jets are restricted to the phase space region where the trans- 
verse momentum of the jet is approximately equal to the photon virtuality, -Px.fwdjet ~ Q 2 ■ This 
condition suppresses the contribution of /^-ordered DGLAP cascades with respect to processes 
unordered in kx B3). The selection of forward jets with a large fraction of the proton energy, 
Xf w djct = Efwdjet/Ep, such that £f w dj C t 3> x, enhances the phase space for BFKL evolution with 
gluon cascades strongly ordered in fractional longitudinal momentum. The above conditions are 
fulfilled by the requirement that the analysed sample contains at least one forward jet which sat- 
isfies the following criteria in the laboratory frame: PT,fwdjct > 6 GeV, 1.73 < ?7f w djct < 2.79, 
a^fwdjct > 0.035 and 0.5 < -P^fwdjet/Q 2 < 6. Here ?7f w djet is the pseudorapidity of the for- 
ward jet. If there is more than one jet fulfilling the above requirements, the jet with the largest 
pseudorapidity is chosen. The upper cut on P| f w djet/Q 2 1S chosen so large in order to reduce 
the contributions of migrations from outside of the analysis phase space, which are due to the 
limited resolution of the Pr.fwdjet measurement. 

The subsample "forward and central jet" is selected by requiring an additional jet in the 
central region of the laboratory frame. This jet is required to have a transverse momentum 
-Pr.cenjet > 4 GeV and to lie in the pseudorapidity region — 1 < 77 cen j e t < 1. The central jet 
must have a large rapidity separation from the most forward jet Ai] = (r/f w dj e t — Vcenjet) > 2. 
This condition enhances the phase space for additional parton emissions between the two jets. 
If there is more than one central jet, the one with the smallest ^ cen j e t is chosen. 

A summary of the selection cuts, defining the DIS phase space for the measurement, the 
forward jet sample and the subsample with an additional central jet, is provided in table[TJ With 
these requirements 13736 and 8871 events are selected for the forward jet and for the forward 
and central jet analysis, respectively. 
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DIS selection 


Forward jets 


Central jets 


0.1 < y < 0.7 
5 < Q 2 < 85 GeV 2 
0.0001 < x < 0.004 


1.73 < ryfwdjet < 2.79 
^T.fwdjot > 6 GeV 
£fwdjct > 0.035 
0-5 < Pl iwd} JQ 2 < 6 


-1 < 7/cenjet < 1 
^T.ccnjct > 4 GeV 
A?7 = 7/f w djet - ^cenjot > 2 



Table 1 : Summary of cuts defining the DIS phase space, the forward jet and the central jet 
selection. If more than one forward jet is found, the jet with the largest ??f w djct is chosen. If there 
is more than one central jet, the one with the smallest ^ ce njet is selected. 

3.3 Cross section determination 



In this measurement in addition to migrations between bins inside the measurement phase space, 
there are considerable migrations from outside of the analysis phase space. This is taken into 
account in the calculation of the cross section corrected to the hadron level: 



a 7 : 



data 



N' 



out 



€i ■ C 



(1) 



Here Nf ata - is the number of observed events in bin i, N° ut is the number of events from outside 
the measurement phase space reconstructed in bin i, and e, is the efficiency in bin i. C is the 
total integrated luminosity. iV° ut and ej are estimated using MC simulations. The purities^] in 
bins of the measured cross sections, as determined from the MC simulations, are at the level of 
80%. 



The efficiency factors e» are calculated according to the formula : 

jydct jyout 



jyhad 



(2) 



where Nf et and A^ had are the numbers of events in bin i at the detector and at the hadron level, 
respectively. For this approach to be valid, the shape of the distributions of all variables on 
which phase space cuts are applied have to be well described by the MC simulations also in 
the phase space extended beyond these cuts. This requirement is found to be satisfied by both 
models considered here. 

The efficiency factors are calculated as the ratio of the model prediction at the detector level 
for a radiative MC and at the hadron level for a non-radiative MC, i.e. the data are also corrected 
for QED radiative effects. The efficiency factors are taken as the average of the factors estimated 
by the RAPGAP and DJANGOH/ ARIADNE models. The uncertainty of the efficiency factors 
is taken to be half of the difference between the factors calculated using the two MC models 
and is included in the systematic error. 



'The purity is defined as the ratio of the number of events generated and reconstructed in the bin to the number 
of events originating from the phase space of the analysis and reconstructed in that bin. 
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3.4 Systematic uncertainties 



The following sources of systematic uncertainties are considered : 

- The model dependence of the bin-by-bin efficiency factors leads to systematic uncer- 
tainties between 2% and 6% for the measured cross sections. 

- The LAr hadronic energy scale uncertainty of 4% for this analysis gives rise to the domi- 
nant uncertainty of 7% to 12% for the measured cross sections. 

- The uncertainty on the electromagnetic energy scale of the SpaCal of 1% results in an 
uncertainty of the measured cross sections below 3%. 

- The uncertainty on the polar angle measurement of the scattered positron of 1 mrad has 
a negligible effect on the cross section measurements. 

- The uncertainty on the determination of the trigger efficiency from the data, using inde- 
pendent trigger samples, leads to an uncertainty between 2% and 4% on the cross section 
measurements. 

- The measurement of the integrated luminosity is accurate to within 1.5%. 

The total systematic uncertainty, adding all individual contributions quadratically, amounts 
to 11 — 12% for the measured cross sections. 



4 Results 

The forward jet cross sections and their uncertainties are given in tableland presented in figures 
(2HU Differential cross sections, da/dA(p, are presented as a function of the azimuthal angle 
difference A<p between the most forward jet and the scattered positron in bins of the variable 
Y = ln(a;f w djet/^)- This variable approximates the rapidity distance between the scattered 
positron and the forward jet. For the selected data sample the normalised shape distributions 
1 jo ■ da/dA(p are also determined, where a is the integrated cross section in a given bin of Y. 
Furthermore, the forward jet cross section is measured as a function of Y. 

The cross section da/dAcf) as a function of A<p is shown in figure[2]for three intervals of the 
variable Y: 2.0 < Y < 3.4, 3.4 < Y < 4.25 and 4.25 <Y< 5.75. These Y bins correspond to 
average x values of 0.0024, 0.0012 and 0.00048, respectively. At higher values of Y the forward 
jet is more decorrelated from the scattered positron. 

The predictions of three QCD-based models with different underlying parton dynamics, 
discussed in section 2, are compared with the data. The cross sections are well described in 
shape and normalisation by CDM which has a BFKL-like approach. Predictions of RAPGAP, 
which implements DGLAP evolution, fall below the data, particularly at large Y. Calculations 
in the CCFM scheme as implemented in CASCADE using the uPDF set AO [El overestimate 
the measured cross section for large A0 values in the two lowest Y intervals. However, this 
model provides as good a description as CDM of the data in the highest Y interval. 
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The shape of the A<p distributions, 1/er • da/dAcf), is compared to the different MC predic- 
tions in the lower part of figure [2l where the ratio R is shown, defined as: 



The precision of the measurements is shown at R — 1 where the statistical and systematic 
uncertainties are indicated. The systematic uncertainty is reduced in the ratio and contains only 
two components added in quadrature: the model dependence of the correction factors and the 
trigger efficiency uncertainty. The ratio plots show that in the analysed phase space region the 
shape of the A0 distributions is well described by all MC models. Since the shape predictions of 
the three models are very similar, this observable alone cannot discriminate among the models. 
It should be noted that the shape of the A<p distributions is rather insensitive to the PDF used 
for event generation. The shape distributions generated using CTEQ6L, CTEQ6M ||3~9l and 
HERAPDF1.0 [|26l differ on average by 1-2%. However, the cross section normalisation is 
more sensitive to the choice of PDF with differences up to 5% for CDM and up to 20% for 
RAPGAP at large Y. 

Predictions of the CCFM model presented in figure [3] indicate a significant sensitivity to the 
choice of the uPDF Set AO and J2003-set 2 give quite different predictions for the differential 
cross sections in all Y intervals. Set AO provides a reasonable description of the measured cross 
sections, except for the region of large A0 in the two lowest Y bins. Predictions using J2003-set 
2 do not describe the data, especially at higher Y, where the estimated cross sections are too 
low. The shape of the A<p distributions is reasonably well described by the set AO. At low Y it 
shows sensitivity to the unintegrated gluon density. 

The cross section da/dY as a function of the rapidity separation Y is shown in figure SI 
The CDM model describes the data well over the whole Y range. The DGLAP predictions fall 
below the data, but approach them at small Y. The predictions of the CCFM model are above 
the data at small Y but describe them well at larger Y corresponding to low values of x. 

The forward and central jet cross sections and their uncertainties are given in table [3] The 
differential cross section da/dA(p as a function of the azimuthal angle difference A0 is shown 
in figure [5] in comparison with the predictions of the three MC models. The cross sections are 
measured in two intervals of Y, 2.0 < Y < 4.0 and 4.0 < Y < 5.75. 

From figure \5\it is observed that at lower Y the predictions of all models describe the cross 
sections reasonably well. At high Y all models undershoot the data: CCFM (set AO) is closest 
to the data, the DGLAP and CDM predictions are below the measured cross section. The ratio 
R in the lower part of figure [5] shows that the shape of the A(p distributions is well described by 
all MC models, as in the case of the forward jet measurements. 

Comparisons of the measured Acp distributions with NLOJET++ predictions are shown in 
figures[6]and[7] The calculations are performed at 0(a 2 s ) precision using the CTEQ6.6 PDF 11301 
and as(M z ) = 0.118. Large theoretical uncertainties of up to 50% from the variation of 
factorisation and renormalisation scales are observed. The size of the theoretical uncertainty 
indicates that in this phase space region higher order contributions are expected to be important. 

In the forward jet sample (figure [6]) for all three ranges of Y the data are above the central 
NLO result but still within the theoretical uncertainty. In the case of the forward and central 




(3) 
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jet sample shown in figure U\ the NLO calculation describes the data at low Y. Only at high Y 
in the regime of the BFKL evolution it is below the data, but again within the large theoretical 
uncertainty. 

In summary, the correlation between the forward jet and the positron decreases with Y and 
the A0 distributions are flat at high Y. The measurements of the forward jet cross sections 
favour CDM and disfavour the RAPGAP model. CASCADE provides a reasonable description 
of the data at large Y, but shows sizeable sensitivity to the uPDF. The shape of the measured 
A0 distributions is well described by MC models based on different QCD evolution schemes. 

The similarity of the A</> shapes of the MC predictions suggests that the forward jet pre- 
dominantly originates from the hard matrix elements which are similar in all three models. 
However, MC studies with RAPGAP show that 80% of the forward jets are produced by parton 
showers. When the initial state parton shower is switched off, the shape of the A</> distribution 
is only slightly changed, but the normalisation is significantly reduced. This indicates that the 
decorrelation in A0 is mainly governed by the phase space requirements, in particular by the 
rapidity separation Y, and that the normalisation of the cross sections is mainly influenced by 
the amount of soft radiation from parton showers, which depends on the evolution scheme. 

5 Conclusions 

Measurements of DIS events at low Q 2 containing a high transverse momentum jet produced 
in the forward direction, at small angles with respect to the proton beam, are presented. Differ- 
ential cross sections and normalised distributions are measured as a function of the azimuthal 
angle difference A0 and the rapidity separation Y between the forward jet and the scattered 
positron. Investigations of the azimuthal correlation between the most forward jet and the out- 
going positron are performed in different regions of Y for the forward jet sample and for the 
subsample with an additional central jet. To test the sensitivity of the measured observables to 
QCD dynamics at low x, the data are compared to QCD models with different parton evolution 
approaches and to predictions of next-to-leading order QCD calculations. 

Measurements of the cross sections as a function of A0 and Y are best described by the 
BFKL-like CDM model, while the DGLAP-based RAPGAP model is substantially below the 
data. The CCFM-based CASCADE provides a reasonable description of the data but shows 
sizeable sensitivity to the unintegrated gluon density. The shape of the A0 distributions does 
not discriminate further between different evolution schemes. The fixed order NLO DGLAP 
predictions are in general below the data, but still in agreement within the large theoretical 
uncertainties. 
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A<p range [rad] 


dcr/dAcj) [pb/rad] 


S s tat [pb/rad] 


Shad [pb/rad] 


S a y S t [pb/rad] 


2.0 < Y < 3.4 


0.0 -0.63 
0.63 - 1.26 
1.26 - 1.89 
1.89 - 2.51 
2.51 - 3.14 


27.3 
33.7 
35.8 
38.9 
47.9 


±3.2 
±3.4 
±3.7 
±3.8 
±4.7 


+2.7 
-2.2 

+2.8 
-3.3 

+3.6 
-4.4 

+4.4 
-4.4 

+4.6 
-3.9 


+1.3 
-1.3 

+1.9 
-1.8 

+2.2 
-1.8 

+2.3 
-2.6 

+2.7 
-2.5 


3.4 < Y < 4.25 


0.0 -0.63 
0.63 - 1.26 
1.26 - 1.89 
1.89 - 2.51 
2.51 - 3.14 


48.2 
56.9 
58.7 
62.9 
60.4 


±4.2 
±4.3 
±4.6 
±4.8 
±4.9 


+5.8 
-4.5 

+6.1 
-6.4 

+7.5 
-6.6 

+6.3 
-6.6 

+7.3 
-7.1 


+2.1 
-2.3 

+2.8 
-2.6 

+2.3 
-2.1 

+2.6 
-3.1 

+2.3 
-2.7 


4.25 < Y < 5.75 


0.0 -0.63 
0.63 - 1.26 
1.26 - 1.89 
1.89 - 2.51 
2.51 - 3.14 


55.1 
60.8 
60.0 
65.0 
57.3 


±4.7 
±5.0 
±4.7 
±5.4 
±5.3 


+6.0 
-5.8 

+5.2 
-6.7 

+7.4 
-7.4 

+7.7 
-7.3 

+5.7 
-4.6 


+2.9 
-3.2 

+2.6 
-2.9 

+4.6 
-4.8 

+4.2 
-4.1 

+4.2 
-4.1 



Y range 


da/dY [pb] 


8 s tat [pb] 


Shad [Pb] 


Ssyst [Pb] 


2.00 - 3.25 


67.9 


±3.3 


+7.5 
-7.7 


+3.0 
-3.1 


3.25-4.00 


194.4 


±6.3 


+21.1 
-20.0 


+8.3 
-8.7 


4.00 - 4.75 


198.2 


±6.7 


+22.5 
-23.4 


+10.4 
-10.4 


4.75 - 5.75 


92.3 


±4.8 


+9.7 
-8.3 


+7.2 
-7.3 



Table 2: Differential forward jet cross section in bins of the variable Y = ln(a;f w djet/^) and the 
azimuthal angle difference A(p between the most forward jet and the scattered positron. The 
statistical uncertainty (S stat ), the uncertainty due to the hadronic energy scale (S had ) and other 
systematic uncertainties (8 syst ) described in the text are given. 
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A<p range [rad] 


dcr/dAcj) [pb/rad] 


Sstat [pb/rad] 


Shad [pb/rad] 


Ssyst [pb/rad] 


2.0 < Y < 4.0 


0.0 -0.63 


18.9 


±2.6 


+1.9 
-1.3 


+1.3 
-1.1 


0.63-1.26 


28.5 


±2.9 


+2.3 
-2.9 


+2.0 
-1.9 


1.26-1.89 


31.6 


±3.4 


+3.9 
-3.8 


+1.9 
-1.9 


1.89 - 2.51 


32.1 


±3.2 


+3.6 
-2.7 


+1.3 
-1.3 


2.51 - 3.14 


33.9 


±3.5 


+2.3 
-3.4 


+2.1 
-2.1 


4.0 < Y < 5.75 


0.0 -0.63 


39.5 


±3.6 


+4.3 
-3.3 


+1.6 
-1.9 


0.63-1.26 


40.8 


±3.6 


+3.4 
-3.9 


+2.2 
-2.2 


1.26 - 1.89 


41.8 


±3.7 


+4.6 
-4.1 


+1.8 
-1.8 


1.89 - 2.51 


43.1 


±4.1 


+5.2 
-4.4 


+2.2 
-2.2 


2.51 - 3.14 


34.9 


±3.7 


+4.0 
-3.7 


+1.7 
-1.7 



Table 3: Differential forward and central jet cross section in bins of of the variable Y = 
ln(xf w djet/^) and the azimuthal angle difference A0 between the most forward jet and the scat- 
tered positron. The statistical uncertainty (5 sta t), the uncertainty due to the hadronic energy 
scale (Sfiad), and other systematic uncertainties (S syst ) described in the text are given. 
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Figure 2: Differential forward jet cross section as a function of the azimuthal angle differ- 
ence A0 between the most forward jet and scattered positron in three intervals of the variable 
Y = m(£f wd j et /£). The inner error bars denote the statistical uncertainties. The systematic 
error due to the uncertainty of the hadronic energy scale is shown separately as a band around 
the data points. Other systematic uncertainties added quadratically to the statistical uncertain- 
ties are represented by the outer error bars. The data are compared with the predictions of 
DJANGOH/ ARIADNE (CDM) and RAPGAP (DGLAP) with HERAPDF1.0, the CASCADE 
predictions (CCFM) are shown with uPDF set AO. In the lower part of the figure the ratio R 
of MC to data for normalised cross sections is shown. The precision of the measurements is 
shown at R = 1 with the statistical and reduced systematic uncertainties indicated as error bars. 
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Figure 3: Differential forward jet cross section as a function of the azimuthal angle difference 
A0 between the most forward jet and the scattered positron in three intervals of the variable 
Y = ln(xf w dj e t/a;)- The data are compared to the predictions of CASCADE (CCFM) with two 
different sets of unintegrated gluon densities. For other details see caption to figure 2. 
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Figure 4: Differential forward jet cross section as a function of the variable Y = ln(xf w dj e t/^)- 
The data are compared with the predictions of DJANGOH/ ARIADNE (CDM) and RAPGAP 
(DGLAP) with HERAPDF1.0, the CASCADE predictions (CCFM) are shown with uPDF set 
AO. For other details see caption to figure 2. 



19 



forward and central jet 



"D 

CO 

B 

& 50 

-e- 
< 

"D 

~6 
"D 



25 



• H1 data 

Energy scale uncert. 

CDM 

CCFM 

DGLAP 



1.25 



0.75 



441 



i 



2.0 < Y < 4.0 
<x> = 0.0018 



< i — — i 



lit! 

i 1 



4 



4.0 < Y < 5.75 
< x > = 0.00058 



R = 



norm. MC 
norm, data 



2 3 

AcKrad) 



Figure 5: Differential forward and central jet cross section as a function of the azimuthal 
angle difference A</> between the most forward jet and the scattered positron in two intervals 
of the variable Y = ln^fwdjct/^)- The data are compared with the predictions of DJAN- 
GOH/ ARIADNE (CDM) and RAPGAP (DGLAP) with HERAPDF1.0, the CASCADE predic- 
tions (CCFM) are shown with uPDF set AO. For other details see caption to figure 2. 
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Figure 6: Differential forward jet cross section as a function of the azimuthal angle difference 
A0 between the most forward jet and the scattered positron in three intervals of the variable 
Y = ln(xf wd j et /a;). The data are compared to the corrected to hadron level NLO predictions 
from NLOJET++ which uses the CTEQ6.6 PDF. Dashed lines above and below the nominal 
NLO prediction show theoretical uncertainty estimated by applying a factor 2 or 1/2 to the 
renormalisation and factorisation scales simultaneously. For other details see caption to figure 
2. 
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Figure 7: Differential forward and central jet cross section as a function of the azimuthal 
angle difference A(p between the most forward jet and the scattered positron in two intervals of 
the variable Y = ln(xf w dj e t/^)- The data are compared to the corrected to hadron level NLO 
predictions from NLOJET++ which uses the CTEQ6.6 PDF. Dashed lines above and below the 
nominal NLO prediction show theoretical uncertainty estimated by applying a factor 2 or 1 /2 
to the renormalisation and factorisation scales simultaneously. For other details see caption to 
figure 2. 
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